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SUMMARY 


A linearized, longitudinal model, in state-space format, is formulated for an 
advanced fighter. The nominal operating point is for a subsonic flight condition. 

The engine is operating with afterburner on. The model is composed of three sub- 
system models: the inlet, the engine, and the airframe. A procedure for combining 

the subsystem models into an integrated model is presented and an integrated con- 
troller is developed by using linear quadratic regulator theory. Notable interaction 
is found in the coupled system. 

A procedure, based on eigenvalue sensitivities, is presented which places the 
feedback gains in a hierarchical arrangement and measures their contribution to the 
optimal solution. With these numbers, called the gain significance matrix, ineffec- 
tual gains can be eliminated thus saving hardware and expense in the realization of 
the physical controller. 


INTRODUCTION 

The desire to produce efficient control systems capable of multivariable input/ 
output and optimal mixing of subsystem interaction has led to programs such as 
INTERACT (integrated research aircraft control technology), AFTI (advanced fighter 
technology integration) , and PROFIT (propulsion- flight control integration technol- 
ogy). Motivation for this work comes from several sources. A primary source is the 
increased interaction between the airplane and propulsion system in certain flight 
conditions. This interaction is intensified with advanced systems and more demanding 
mission requirements. Another incentive is the need to improve efficiencies. As 
component efficiencies become more difficult to improve, the alternative is to 
iitprove overall efficiency through optimal integration and control of conponents. 
Finally, the historical separation of major airplane systems in the design process 
has not fully exploited the capabilities of the overall system; this is iirportant in 
light of the wide operational envelope of a modern -fighter. 

The benefits of integrated control design have been demonstrated. For example, 
Michael and Farrar (ref. 1) developed an integrated inlet-engine controller for the 
Pratt and Whitney F401 engine and an internal conpression supersonic inlet. They 
were able to improve steady- state thrust by 6 percent and reduce the variation of 
normal-shock position by a factor of three. This predicts a substantial iitprovement 
in performance of the inlet- engine system. 

In the same spirit, this report contributes to integrated design methodology by 
achieving two objectives. First, a state-space model of an airplane comprised of 
three major subsystem models (inlet, engine, and airframe) is provided as a working 
tool. A procedure is shown for combining the subsystem models into an integrated 
model, and then an integrated controller is developed by using linear-quadratic regu- 
lator (LQR) theory. The second objective is to provide a method to discern the rela- 
tive importance of the feedback gains by conputing a gain significance matrix, which 
is a function of the eigenvalue sensitivity to the gains. This matrix indicates the 
ineffectual gains which can be eliminated; thus, hardware and expense may be saved in 



the realization of the physical controller. It also indicates sensitive gains so 
that beneficial interactions can be identified and adverse interactions can be 
attenuated. 


SYMBOLS 
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w 


B 
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ik 
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D,I 


D 


D 

F 


F,G 


kJl 

H 

h 

I 

J 

^a2 

k f i f j 


uncoupled plant matrix 

coupled plant matrix 

coupled, closed-loop plant matrix 

2 

xnlet capture area, ft 

, 2 
jet area, ft 

2 

wing area, ft 

control distribution matrix for uncoupled system 
control distribution matrix for coupled system 
element of B matrix 

output equation state matrix for uncoupled system 
output equation state matrix for coupled system 
drag coefficient for inlet 
moment coefficient for inlet 

output equation control matrix for uncoupled system 

output equation control matrix for coupled system 

feedback gain matrix 

coupling equation matrices 

element of feedback gain matrix 

output equation matrix 

altitude, ft 

unit (identity) matrix 

performance index 

distortion factor 

matrix and vector elements 
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LQR 

M 

m 

'"c 

mf 

"c 

PLA 

Pst 

Pt 

^t,0 

Pt,2 

Pt,2/Pt 


Q 

q 

R 


1/s 

T 


T 


f 


Tfrj, 

t 


^1/2 

U 


U 


U 


n 


linear-quadratic regulator 
Mach number 

number of states or eigenvalues and eigenvectors 

compressor surge margin, percent 

fan surge margin, percent 

compressor speed, rpm 

fan speed, rpm 

power lever angle, deg 

static pressure, psi 

total pressure, psi 

nominal total pressure, psi 

engine face total pressure, psi 

pressure recovery, ratio of engine face total pressure to free- stream total 
pressure 

output weighting matrix in performance index 
pitch rate, rad/sec 

control weighting matrix in performance index 
gain significance matrix for the nth eigenvalue 
integrator 
thrust, lb 

fan inlet temperature, 

fan turbine inlet tenperature, °R 

time, sec 

time to damp to one-half amplitude, sec 
time to daitp to double amplitude, sec 
control input vector for uncoupled system 
control input vector for coupled system 
nth right eigenvector for closed-loop system 
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V 


velocity, ft/sec 


nth left eigenvector for closed-loop system 
^a,E engine airflow, Ib/sec 

(^a,E> engine airflow trim request, Ib/sec 

w ^ inlet airflow, Ib/sec 

fuel flow, Ib/hr 

^f,A/B afterburner fuel flow, Ib/hr 
X state vector 

X state vector differentiated with respect to time 

throat ramp position 

Y output vector 

a angle of attack, rad 

P parameter in gain significance analysis 

\ 

Y flight-path angle, rad 

6^ elevator control, rad 

e 

C damping coefficient 

0 pitch attitude^ rad 

X eigenvalue 

inlet duct time constant, sec 
throat ramp time constant, sec 
rotating ramp time constant, sec 
<l>^ rotating ramp position, deg 

to natural frequency, rad/sec 

Subscripts: 

A airframe 

C commanded 

E engine 

1 inlet 
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n 


nth 


sp short period 

Superscript: 

T transpose 


MODEL DEVELOPMENT 

The linear time-invariant airplane model developed in this report is synthesized 
from three linear subsystem models. The subsystem models are the inlet, the engine, 
and the airframe aerodynamic models. The integrated system is a model of an advanced 
fighter with twin turbofan engines. The flight condition about which the models are 
linearized is subsonic and afterburners are in use. Table 1 summarizes the nominal 
operating values for most state and output variables. 

The airframe model contains only longitudinal aerodynamics. Despite the lack of 
lateral dynamics, there is still notable interaction among the three subsystems as 
will be explained in the section "Discussion.” The aerodynamic data for the airframe 
model were generated with an operating-point program based on the data in refer- 
ence 2. The airframe model, in state-space format, is given in table 2. 

The engine and inlet models were taken from reference 3. The engine and inlet 
models, in state-space format, are given in tables 3 and 4. The linear inlet model 
is a two-dimensional mixed compression inlet. Since the inlet model was substan- 
tially modified from that given in reference 3, a block diagram of the present* inlet 
model is provided in figure 1 to allow a direct comparison with the original model in 
reference 3. The primary modification to the inlet model is the addition of j 
and T^ as output variables. The inlet moment coefficient j was estimated with 
the operating-point program mentioned in the previous paragraph, and the fan inlet 
temperature T^ was based on the model in reference 4. 

The synthesis of the component models is readily accomplished after each is put 
into state-space format. The analyst needs only to specify the coupling equation 
which completely defines the subsystem interactions. The coupling equation is 


U 


GU + FY 


(1) 


where U is the control input vector for all three subsystems, Y is the output 
vector for all three subsystems, F and G are the coupling equation matrices, and 
U is the input vector for the integrated (coupled) system. The U and Y vectors 
have the form 


< 

II 

D 


«e]" 

(2a) 

Y = [^A 



(2b) 
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The vector U contains the system inputs which are independent of any system output. 
Thus, equation (1) is simply a mathematical statement of the input/output relation- 
ships among the subsystems, Ihis relationship is shown diagrammatically in fig- 
ure 2, The G and F matrices are given in table 5, 

The integrated system is obtained by substituting the coupling equation 
(eq, (1)) into the following state-space equations for the three subsystems: 


X = AX + BU (3a) 

Y = CX + DU (3b) 

In equations (3), U and Y are given by equations (2) and X is in the same for- 
mat as equations (2): 

X = [^A (4) 

The matrices A, B, C, and D correspond in the appropriate way; for example, the 
block diagonal matrix A is 


A = diag[A^ A^ A^] 


After 

by 

substitution of equation (1) into equations (3), 

the integrated system is given 


, A /N 

X = AX + BU 

(6a) 


A A A 

Y = CX + DU 

(6b) 

where 

A 

X, Y, and U are as indicated before and 



A = A + BF(I - DF)“^C 

(7a) 


B = BG + BP(I - DF)“^DG 

(7b) 


C = (I - DF)~''c 

(7c) 


D = (I - DF)“''dG 

(7d) 


Matrices A, B, C, and D are also presented in table 6 . 
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The solution to equations (6) is accomplished by using standard LQR theory. In 
this report^ ORACLS software was used from reference 5. The ORACLS library provides 
routines to solve the regulator problem of the form 


X = AX + BU 


(8a) 


Y = HX (8b) 

J = / (X*^QX + u"^RU) dt (8c) 


Since the desired problem is not in this form, another control variable transforma- 
tion is required to eliminate the cross-product terms (eq. (10b)) in the performance 
index. This procedure is covered in reference 5. The desired form is 


X = AX + BU 
Y = CT + DU 
J = /(y'^qy + u'^ru) dt 


(9a) 

(9b) 

(9c) 


Equation (9c) can be expanded to show the cross-product terms by substituting equa- 
tion (9b) into equation (9c) to get 


J = /[(ct + DU)‘^Q(CX + Du) + U'^RU] dt 


(10a) 


and expanding gives 


J = /[x'^c'^QCX + 


2X^C*^QDU + U^(r + D-^QDJUJ dt 


(10b) 


Ihe final statement of the LQR problem is the specification of the weighting 
matrices. These matrices are basically design parameters chosen by the control 
designer to achieve the desired response characteristics. The criteria for choosing 
these matrices are primarily based upon experience. Since methodology is more impor- 
tant than strict model fidelity for this report, two simple criteria were used to 
select weights. The first, suggested in reference 6, takes the ith weight to be the 
inverse of the squared maximum allowable deviation of the ith parameter. The second 
required that the short-period frequency and damping meet the Military Specification 
requirements (ref. 7) for the category and class of airplane. The first criterion 
was used to obtain initial estimates of the weighting matrices. The second criterion 
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was satisfied through a trial-and-error process of adjusting the initial weights. 
The final weighting matrices Q and R are given by the following equations: 


Q = diag[1.0E-05 0. 0.1 0.25 0. 0.4 9. 0.01 1. 1. 0. 0.03 0.04 0.04 0.002] 


R = diag[4. 0.04 0.02 1.] 


FEEDBACK GAIN ANALYSIS 

The solution to the LQR problem for the integrated system is a feedback gain 
matrix. This matrix defines the optimal control law with full state feedback. It is 
desirable to eliminate ineffectual gains since they represent added weight and hard- 
ware in the physical system. Of course, eliminating any element of the gain matrix 
means the solution is no longer the optimal solution; however, the solution may still 
be close enough to optimal from an engineering point of view. With this in mind, a 
method is shown that will determine which gains are essential, which gains can be 
eliminated, and in what order they should be eliminated. In addition, a relative 
measure of the penalty for removing the gain is provided. 

The procedure is based upon the closed-loop eigenvalue sensitivity to gain 
matrix elements. Sensitivity problems have been considered by many authors. (For 
example, see refs. 8 through 12.) The sensitivity of the nth eigenvalue to the 

gain matrix element f^j^ is given by 


ax 


kl 


m 


i±=^ 




u 


n, X 


( 11 ) 


where and are the^nth left and right eigenvectors of the closed-loop- system 

matrix, respectively, and b^]^ is the (i,k) element of the control input distribu- 
tion matrix. This result is developed in the appendix. 

In order to compare the sensitivities in a hierarchical arrangement, a gain sig- 
nificance matrix is defined for each eigenvalue. For the nth eigenvalue, the 

(k, JO element of is given by 


^n,kJl 


ax 


n 


5f. 


kJl 



( 12 ) 


The elements of this matrix represent the nondimensional modulus of the sensitivity 
of the nth eigenvalue to each gain matrix element. This can be interpreted as a 
relative indication of the penalty for changing any particular gain and a hierarchy 
for eliminating gains. In general, if for all eigenvalues is much less 

than 1, the removal of the corresponding gain f^^^ will not appreciably affect the 
optimal solution. If for any eigenvalues is greater than 1, the correspond- 

ing gain element should be considered essential. 
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DISCUSSION 


The integration of the three major airplane subsystems (airframe, inlet, and 
engine) primarily affected the airframe stability characteristics. Ihe effects on 
the inlet and engine modes were negligible. Table 7 shows the open-loop eigenvalues 
for the coupled and uncoupled systems and each mode is identified with its corre- 
sponding subsystem. The aerodynamic model is fifth order; the five states are V, 
a, q, 0, and h. As shown in table 7, the oscillatory short-period mode, involving 
primarily a and q, has changed from o^n^sp ^sp == ^,sp == 

Csp = 0.39 after coupling the subsystems. Ihe phugoia mode, involving primarily 
V and 0, becomes oscillatory and remains unstable after integrating the subsystems. 
The phugoid natural frequency is 0.054 rad/sec and the time to damp to double ampli- 
tude t 2 is 23.6 sec. The last mode of the airplane model is associated with the 
altitude h. It has a stable response both before and after the systems are com- 
bined. The times to damp to one-half amplitude before and after coupling the 

systems are 10.1 and 7.2 sec, respectively. 

After integrating the subsystems with coupling equation (1), the standard LQR 
problem was solved. The open- and closed-loop eigenvalues along with an indication 
of what subsystem they correspond to are given in table 8. The results show that the 
inlet modes were virtually unaffected by the addition of feedback. The engine modes 
were slightly changed primarily through increases in damping, although the mode asso- 
ciated with fan speed had a more substantial change. This mode changed its time to 
damp to one-half amplitude from t ^^2 = sec (open loop) to t^ ^2 = sec 

(closed loop). The aerodynamic model was the most affected by feedback. The short 
period mode changed its natural frequency from sp ^,sp == 

damping ratio changed from Csp = 0*39 (open loop5 to Csp = ^(closed loop). 

Of course, this was designed into the model by the appropriate choice of weighting 
matrices Q and R in the performance index. The phugoid mode became stable with 
feedback and the eigenvalue associated with altitude changed from -0.096 (open loop) 
to -0.053 (closed loop). 

The resulting airplane model has enough interaction between the propulsion sub- 
systems and the airframe subsystem to demand cross feedback in the solution to the 
LQR problem. This provides a useful tool for testing the gain significance matrix 
as a measure of feedback gain importance to the, optimal solution. The values 
of for this model are given in table 9. To demonstrate the utility of the gain 

significance matrix, gain elements were eliminated according to the following rule: 
if the (k,D element of all the is less than p, remove the (k,X) element of 

the feedback gain matrix. The p parameter was chosen to be 0., 0.01, 0.1, 0.5, 1., 
and 00 . Choosing p = 0. results in no elements of the gain matrix being elimi- 
nated. This is the full state feedback case and the optimal solution to the LQR 
problem. Choosing p = oo eliminates all the gains and so represents the open-loop 
solution. The closed-loop eigenvalues for each value of p are given in table 10. 
The location of the short-period pole for varying p is given in figure 3. As men- 
tioned before, when an element of is close to unity, the corresponding gain 

element is essential for the optimal solution. This is demonstrated in table 10 and 
figure 3 since negligible change in the eigenvalues occurs for p < 0.01, even though 
25 percent of the gains have been eliminated. In addition, with p < 0.1, 65 percent 
of the gains are eliminated while still maintaining acceptable performance. For 
larger values of p, the change in eigenvalues becomes much larger and the system 
designer must choose the acceptable trade-off between closeness to the optimal solu- 
tion and reduction in the nvimber of feedback gains. In addition, for large values of 
P, the possibility exists for eliminating certain combinations of gains which may 
cause an air^lif ied adverse response; this is shown in figure 3 for p = 1 . 
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The elements also provide a relative measure of the penalty for removing a 

particular gain. It is only a relative measure since eigenvalues and their sensi- 
tivities to gains can be coitplex numbers. In simplifying these numbers with the 

formula/ only the modulus is used; therefore, phase information is lost (relative 
size coitparisons make sense only on the real axis) . Therefore, the elements 

cannot be used directly to compute the actual change in an eigenvalue after modifying 
a gain. This is a penalty for simplification. The can still be a useful indi- 

cator of the cost of changing or removing a gain element. As was seen in table 10(b) 
for p = 0.01, very small elements indicate very small or negligible changes in 

the associated eigenvalue when the corresponding gains are removed. Now, for exam- 
ple, consider element (1,3) for eigenvalues 7 and 8, the short-period mode 

(table 9); this gives , 3 ) » 0.69. This corresponds to the gain element which 

feeds back pitch rate q to the elevator control 6^. The result of removing only 
this gain should be significant. The resulting closed-loop eigenvalues are given in 
table 11. The change in natural frequency is from w^^sp “ ^,sp ~ 5.6, and 

the change in damping ratio is from C = 0.64 to C - 0.13, almost a factor of 5. 
These changes show that indicates where the sensitive and more important gain 

elements are for each mode, although it does not indicate whether the sensitivity 
will be expressed in the magnitude or phase of the eigenvalue. 


CONCLUDING REMARKS 

This report satisfies two objectives. The first objective was to develop a 
simple model of a modern fighter with interactive propulsion and airframe subsystems. 
This provides a working tool for studying integrated control methodology. The second 
objective was to develop a method for eliminating ineffectual feedback gains. This 
provides a method for simplifying the control system without major performance 
penalties. 

There are limitations to these two efforts, however. Ihe airplane model has 
only longitudinal aerodynamics and is linearized about one flight condition. The 
engine is operating with afterburner on. This flight condition represents only one 
small area in the fighter's flight envelope. The gain significance matrix pro- 

vides an easily implemented measure of the relative importance of the feedback gains 
and an indirect measure of the actual change that can occur in the eigenvalues. The 
actual change in the eigenvalue can be calculated from the eigenvalue sensitivities 
if these are retained during the calculation of 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
July 15, 1982 
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APPENDIX 


DERIVATION OF EIGENVALUE SENSITIVITIES TO FEEDBACK GAINS 

The derivation of eigenvalue sensitivities to feedback gains closely follows the 
development for sensitivities to the system matrix, (See ref, 8 ,) 

Consider the mth order system matrix A, the closed-loop system A with dis- 
tinct eigenvalues, the distribution matrix B, and the gain matrix F where 


A = A + BF 


(A1) 


Let Vj and be left and right eigenvectors of A so that 


“i = 


(A2) 


and 




(A3) 


and 


m m 

ufv. = v:u, = 6. . 

ID D 1 3-D 


(A4) 


where 6 is the Kronecker delta. 


Differentiating equation (A2) with respect to gain 


~ dU. 5 \a dU. 

U 4 + A y j ;' " ~ — U' + 


■^kA 


kA 


(A5) 


Multiplying equation (A5) on the left by vt gives 


„T 9A „ , „T~ ^^i _ „T ^^i „ , „T, ^“i 


kA 


(A6) 


Substituting equation (A3) into equation (A 6 ) and dropping the i subscript gives 


yT 9A _ u 
^^kA 



(A7) 
1 1 



APPENDIX 


Note that 



(A8) 


Substitute equation (A8) into equation (A7) to get 


kA 


v'^B 


5F 


W 


u 


kJl 


(A9) 


or 



vTS6ij,6j^U = 


(A10) 


Thus the kth element of the vector v"^B times the Ath element of vector U is the 
eigenvalue sensitivity to the gain Equation (A10) can also be written as 


df^ ” ( 2 ^n,i^ik)^n,A 


(All) 


which emphasizes that the equation is for the nth eigenvalue. 
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Airframe 


Inlet 



Engine 


n^ = 10 040 rpm 
n^ = 12 930 rpm 
= 8182 Ib/hr 
'^f,A/B = 24 870 Ib/hr 
T = 16 930 Ib/engine 
w, „ = 173.2 Ib/sec 
PLA = 120® 

Aj = 4.7 ffc2 

= 19.35 percent 
= 20.39 percent 
Tf^ = 2171 »r 
( w^^E^tr ~ 225 Ib/sec 










m fO fH O O O rH CV O O 0^0 000 


TABLE 2.- AIRFRAME MODEL 


Xj . [v 

Xft = [m a q Y 
“a " 2’' 2C„_j 2C„_jP 


a q 0 h]^ 

q Y h]T 


^A = Va + ^A^A 
""a = C^^A °A^A 


matrix 

•397P-02 

.<J15E-C'5 

.324E-03 


matrix 

•223E+00 

.910E-M 

•483E+01 


matrix 

•646E-04 


matrix 


•4.037E + 01 
•1.5C>5E+00 
•1.292E+01 

0 • 

•9.332E + 02 


P.9!58E-04 

•?.ei2E-O0 

1.261E-0B 


0 . 

0 , 


0 • 

l.OfeCE +00 

0 . 

•l.OOCE+CO 

0 . 


D. =lo( 


0 • 

l.OOJE+tO 

•2.142E+00 

1 . 0 i>aE+t 0 

C. 


•3.005E+02 

C. 

0 , 

0 * 

c. 


0 , 

0 . 

1 .OOOE+00 

c . 

0 • 


•3.217E+01 

0 . 

0 . 

0 . 

9*332E+02 


0 « 

0 . 

3.207E+01 

0 . 

0 . 


0 * 

0 . 

0 . 

l.UOOE+00 

0 . 


1.32PE-04 

1.493E-06 

7.503E-06 

C* 

c. 


c« 

0 . 

0 . 

0 • 

I .ODCE-fbd 



TABLE 3.- ENGINE MODEL 


^E = '"a,E "‘c ’^fT.^'^ 

Ug=[PLA (w^,E^tr Pt,2 

^E = ^^E + ®eUe 


Ag matrix 

■3.385E-01 -1.0P9E-01 

3.318E-01 -1.270E+00 

•l,462E+tO -<»,795E+00 

•1.362E + 01 6.116E+01 

Bg matrix 

1.1«8E+01 6.292E+01 

2.254P+00 6.3B9E+00 

-V.^P6E + 00 -3.679E + C1 

1.118E+04 -1.222E+02 


Cg matrix 

-1.281E+00 

2.029E-02 

2.485E-02 

-1.595E-03 

9.524E-02 


9.396E-01 

-1.167E-02 

-1.756E-02 

5.036E-03 

-2.514E-02 


= ^^E ^ °E«E 


1.269E-01 

3.654E-01 

-4.503E-01 

1.214E+01 


7.687E+01 

2,575E+02 

1.067E+03 

3.526E+03 


1.147E+00 

1.045E-03 

-7.241E-03 

-8,618E-04 

7.329E-02 


8.939E+01 

1.371E+01 

7.705E+00 

1.091E+00 

7.130E+01 


-2.394E-03 

-3.013E-03 

-3.114E-03 

-7.291E+00 


1.330E>00 

1.462E+01 

4.153E+01 

1«043E<^02 


1.502E-01 

9.343E-05 

6.049E-05 

1.122E-05 

-1.879E-04 


-5.816E+01 

-2.416E-01 

-2.665E-02 

-6.044E-02 

2.149E+00 


Dg matrix 

2.036E+01 -2.690E+01 
5.862E-02 7.717E-02 
1.643E-02 3.773E-C1 
1.779E-02 -2.870E-C2 
8.597E-01 1.371E+00 


Pst '',2 ^j]’' 


0. 

0. 

0. 

0. 


0. 

0. 

0 . 

0. 


“8.769E+02 

0. 

0. 

0. 

0. 


0 . 

0. 

1.030E+01 

0 . 

0. 


3.256E+02 

5.239E+01 

-1.214E+03 

-1.440E+03 


-1.564E+03 

0.646E+OO 

1.863E+01 

-7.235E-02 

B.600E<^00 



moo o o o 


TABLE 4.- INLET MODEL 



Aj matrix 

•125E+01 0. 

. -3.125E+01 

• 0 . 

Bj matrix 

. 1.791F+03 

. 2.472E+01 

• Of 

Cj matrix 

-6.8C0E-03 5.475P-01 

2*200E-02 -3,<»90E+00 

0 • 0 • 

-3.500E-05 -o.lOOE-0^ 
-1.063E-03 0. 

Dj matrix 

1.215E + 01 -<».583E-01 

2.0A0E-C1 1.500E+00 

6.320E+01 0. 

1.150E-03 2.570E-03 

7.300E-03 -8.600E-03 


0 . 

0 . 

-2.857E+01 


0 . 0 . 

0 . 0 . 

0, 2.857E+01 


3.430E-03 

8.000E-03 

0 . 

-1.400E-05 

-4.660E-05 


-4.050E-04 0. 

0 . 0 . 

-1.^705-02 0. 

0 . 0 . 
0 . 0 . 
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TABLE 5.- G AND F MATRICES FOR COUPLING EQUATION 
[Matrix order is indicated by first row elements] 


Coupling equation: U = GU + FY 
G matrix 


ItOOOE^OO 

0. 

0 . 

0. 

0. 

0. 

0* 

0. 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 


0. 

0. 

0« 

c • 

0. 

0, 

0. 

0* 

1. 000^^00 

0. 

0. 

0. 

0. 

0 . 

0* 


0. 

0. 

0 . 

0 . 

0. 

0. 

c. 

0* 

0. 

I.OOOE-t-CO 

0. 

0 . 

0« 

0 . 

c. 


0 . 

0. 

0. 

c. 

0. 

0. 

0 . 

0« 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

l.OOOE^OO 


A 

F matrix 


( 1 , 2 ) 

0. 

0 . 

0 . 

0« 

0 . 

l.OOOE+00 

0 . 

0. 

0. 

0, 

0 . 

0 . 

0. 

0. 

0* 


(1, 5) 

0. 

0. 

0. 

0. 

0* 

c • 

X.OO0E400 

0. 

0. 

0. 

0. 

0. 

-2.830E-0^ 

0* 

0. 


( 1 , 6 ) 

0. 

0. 

0* 

0. 

0. 

0. 

0. 

0. 

0. 

0* 

1.000E400 

0* 

0. 

0. 

0 . 


(1, 7) 

0. 

0* 

0 . 

0. 

0. 

0 . 

0. 

0. 

0, 

0. 

0 . 

0 . 

0 . 

1«OOOE^OO 

0 . 


0. (1’^) 

0 • 

2«000E^00 

0* 


0. <*•“' 

t. 

0 . 

2»OOOE^OO 

0. 

0. 

0 . 

0 . 

0. 

0. 

0 . 

0 . 

0. 

0« 

0, 


0 . <*•“> 

2.000E+00 

0. 


0 . 

0 . 

0 . 

0 . 

0. 

0« 

0. 

l.OOOE-^00 

0 . 

0. 

0. 

0. 

0. 

0 . 

0 . 


(1, 13) 


(1, 14) 
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TABLE 6.- INTEGRATED MODEL 


Matrix order is similar to table 5 
X = [V a q 0 h 

Y = [m a q Y Pt,2 ^a2 "^fT ‘^D,! *^M,I 

G=[6^ PLA (w^^E^r ^j]’’ 

X = AX + BU 

Y = CX + DU 


A matrix 


-4.112E-02 

-A.199E+01 

0. 

-3.217E+01 

2.253E-03 

1.994E-02 

6.34?5E-01 

-■9.895E-05 

-1.505E+00 

l,000E+00 

0. 

1.427E-06 

3.418E-08 

-2.7f5?E-06 

-8.727E-04 

-1.347E+01 

-2.142E+C0 

C. 

7.533E-06 

-6.818E-02 

1.234E-06 

0. 

0, 

l.OOOE+00 

0, 

C« 

0. 

0. 

0. 

-9.332E+02 

0. 

P.332E+02 

0. 

0. 

0. 

3.858E-14 

1.791E+03 

0. 

0. 

-1.295E-15 

-3.125E+01 

2.407E-11 

5.326E-16 

2.472E+01 

0. 

0, 

-1.78FE-17 

-2.252E-15 

-3.125E+01 

A.037E+00 

-1.795E+02 

C. 

c. 

-4.336E-02 

-2.664E+00 

2.148E+02 

1.008E+00 

-3.^23E+01 

C. 

0. 

-5.334E-02 

-5.227E-01 

4.209E+01 

-1.8^4E+00 

1 .180E+02 

0. 

0* 

-1.399E-C1 

1.751E+00 

-1.410E+02 

9.172E+00 

-4.890E+02 

0. 

0. 

2.614E-01 

-7.256E+00 

5.842E+02 

A,969E+0i 

-1.616E+03 

0* 

0. 

-3.170E+00 

-2.398E+01 

1.930E+03 

8.962E-03 

-2.295F-C3 

1.683E-C3 

2.055E-03 

2.691E-04 



-l,724E-08 

7,203F-0P 

-5.28^E-08 

-6.45CE-08 

“8.446E-09 



-2.989E-03 

-3.231E-08 

2.370E-08 

2.893E-08 

3.788E-09 



0. 

0. 

0. 

0. 

0. 



0. 

U • 

0. 

U* 

0. 



-3.921E-14 

c. 

C, 

0. 

0. 



-5.413E-16 

0. 

0. 

0. 

0. 



-2.723E+01 

5.797E-C1 

-3.334E-C1 

2.986E-02 

2.669E-03 



2.637E-01 

-3.385E-01 

-1.089E-01 

1.269E-01 

-2.394E-03 



-8.832E-01 

3.318E-C1 

-1.270E+00 

3.654E-01 

-3.013E-03 



3.660E+00 

-1.462E+C0 

-4.795E+00 

-4.503E-01 

-3.114E-03 



1.209E+01 

-1.362E+C1 

6.116E401 

1.214E+01 

-7.291E+00 






to 

o 

A 

B matrix 

TABLE 6.- 

Continued 



-9.223E+00 

-3.646E-02 

-4.819E-C2 

-2.802E+00 


-1.910E-01 

1.145E-06 

l,513E-06 

8.795E-05 


-2.483E+01 

-5.13fE-07 

-6,784E-C7 

-3.944E-05 


0. 

0* 

0. 

0. 


0. 

0. 

0. 

C. 


0. 

0. 

0. 

0. 


0. 

0 • 

t< • 

0. 


0. 

-1.67?5E + 00 

2#205E+C0 

2.470E+02 


0. 

1.148E+01 

6.292E+D1 

3.256E+02 


0. 

2,254E+00 

8.389E+10 

5.239E+01 


0. 

-4.486E+00 

-3.679E+01 

-1.214E+03 


0. 

1.118E+04 

-1.222E+C2 

-1.440E+03 


A 




D matrix 




0. 

0. 

0, 

0. 

0. 

0, 

0. 

C. 

0. 

0, 

0« 

0. 

C« 

o* 

G* 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

0 • 

c. 

0. 

0. 

-2.036E+01 

~2.690E+01 

"1.564E+03 

0. 

-5.862E-02 

7.717E-02 

8.64EE+00 

0. 

1.643E-02 

3.773E-C1 

1.863E+01 

0. 

-1.779E-02 

-2.870E-02 

-7.235E-02 

0* 

8.597E-01 

1.371E+CO 

8.800E+00 



TABLE 6#- Concluded 


/V 

C matrix 


9.646E-04 

2.155E-17 

0. 

0 . 

0. 

1.172E-Il)2 

1.968E-04 

0.O25E-O2 

1.109E-06 

7.041E-06 

-3.620E+00 

1.413E-01 

8.613E-02 

7.936E-03 

-6.631E-01 


3,808E-19 

-2.190E-17 

0. 

0, 

0. 

3.430E-03 

8.000E-03 

2.437E-17 

-1.400E-05 

-4.660E-05 

3.066E-01 

4.703E-02 

-5.597E-02 

3.742E-03 

-2.446E-01 


-5.551E-17 

l.OOOE+OO 

0 • 

-l.OOOE+OO 

0 . 

-4.5P3E-01 

1.500F+JC 

-3.553E-15 

2.570E-C3 

-8.600E-03 

-4.097E+01 

-6.2P3E+00 

-1.898E+01 

-5.0C0E-01 

3.268E+01 


0. 

0 . 

0 . 

C. 

0. 

c. 

0 . 

0 . 

0. 

-1.281E+C0 

2.029E-C2 

2.485E-C2 

-1.595E-C3 

9.524E-02 


0 . 

0 1 

l.COOE+OO 

0 . 

0. 

0 . 

c. 

0 . 

c • 

0 . 

c • 

0 . 

Cl. 

0 . 

0 . 


0. 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 

0 . 

0 . 

0 . 

9.396E-01 

-1.167E-02 

-1.756E-C2 

5.036E-03 

-2.514E-C2 


0 . 

0 . 

C. 

l.OOOE+OO 

0. 

0 . 

c. 

0. 

0 . 

0 . 

0 . 

c. 

0. 

0 . 

0. 


0 . 

0 . 

0. 
c. 
c . 

0. 

0 . 

0 . 

0 . 

0 . 

1.147E+00 

1.045E-03 

-7.241E-03 

-e.618E-04 

7.329E-C2 


1.735E-18 

-7.233E-19 

0 . 

0 . 

l.COOE+CO 

-4.050E-04 

-2.168E-10 

-1.670E-02 

-5.082E-21 

2.711E-20 

1.183E+00 

-1.518E-03 

-2.675E-03 

5.675E-04 

-7.012E-C3 


0. 

0 . 

0. 

0 . 

0. 

0. 

0 . 

0 . 

0 . 

0 . 

1.502E-01 

O.343E-05 

6.049E-05 

1.122E-05 

-1.879E-04 


-7.55CE-19 

-9.109E-17 

0 . 

0. 

0 • 

-6.800E-03 

2.200E-02 

-4.832E-17 

-3.500E-05 

-1.063E-03 

-6.079E-01 

-9.323E-02 

-2.790E-01 

-7.419E-03 

4.848E-01 


6.078E-17 

1.344C-14 

0 . 

0. 

0 . 

5.475E-01 

3.490E+00 

3.89CE-18 

9.100E-04 

9.023E-1R 

4.894E+01 

7.506E+00 

4.017F+01 

5.973E-01 

3.904E+01 



TABLE 7.- OPEN-LOOP EIGENVALUES FOR COUPLED AND UNCOUPLED SYSTEMS 


Subsystem 

Re(\) 

lm( X) 

0) 

c 

ti /2 or +2 

Open-loop eigenvalues for uncoupled system 

Airframe 

9.350 E-03 

0. 



7.41 E+01 

Airframe 

3.109 E-02 

0. 



2.23 E+01 

Airframe 

-6.895 E-02 

0. 



1.01 E+01 

Engine 

-4.938 E-01 

0. 



1.40 E+00 

Engine 

-7.985 E-01 

1.313 E+00 

1.537 E+00 

5.195 E-01 

8.68 E-01 

Engine 

-7.985 E-01 

-1.313 E+00 

1.537 E+00 

5.195 E-01 

8.68 E-01 

Airframe 

-1.826 E+00 

3.582 E+00 

4.021 E+00 

4.542 E-01 

3.79 E-01 

Airframe 

-1.826 E+00 

-3.582 E+00 

4.021 E+00 

4.542 E-01 

3.79 E-01 

Engine 

-7.259 E+00 

0. 



9.55 E-02 

Inlet 

-2.857 E+01 

0. 



2.43 E-02 

Inlet 

-3.125 E+01 

0. 



2.22 E-02 

Inlet 

-3.125 E+01 

0. 



2.22 E-02 

Open-loop eigenvalues for coupled system 

Airframe 

2.941 E-02 

4.471 E-02 

5.352 E-02 

5.496 E-01 

2.36 E+01 

Airframe 

2.941 E-02 

-4.471 E-02 

5.352 E-02 

5.496 E-01 

2.36 E+01 

Airframe 

-9.664 E-02 

0. 



7.17 E+00 

Engine 

-4.443 E-01 

0. 



1.56 E+00 

Engine 

-8.106 E-01 

1.314 E+00 

1.544 E+00 

5.250 E-01 

8.55 E-01 

Engine 

-8.106 E-01 

-1.314 E+00 

1.544 E+00 

5.250 E-01 

8.55 E-01 

Airframe 

-1.757 E+00 

4.174 E+00 

4.529 E+00 

3.880 E-01 

3.95 E-01 

Airf raime 

-1.757 E+00 

-4.174 E+00 

4.529 E+00 

3.880 E-01 

3.95 E-01 

Engine 

-7.259 E+00 

0. 



9.55 E-02 

Inlet 

-2.725 E+01 

0. 



2.54 E-02 

Inlet 

-3.125 E+01 

0. 



2.22 E-02 

Inlet 

-3.139 E+01 

0. 



2.21 E-02 
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TABLE 8.- OPEN- AND CLOSED-LOOP EIGENVALUES FOR COUPLED SYSTEM 


Subsystem 

Re(\) 

lm( X) 

03 

C 

^1/2 ^2 

Open- loop eigenvalues for coupled system 

Airframe 

2.941 E-02 

4.471 E-02 

5.352 E-02 

5.496 E-01 

2.36 E+01 

Airframe 

2.941 E-02 

-4.471 E-02 

5.352 E-02 

5.496 E-01 

2.36 E+01 

Airframe 

-9.664 E-02 

0. 



7.17 B+00 

Engine 

-4.443 E-01 

0. 



1 .56 E+00 

Engine 

-8.106 E-01 

1.314 E+00 

1.544 B+00 

5.250 E-01 

8.55 E-01 

Engine 

-8.106 E-01 

-1.314 E+00 

1.544 E+00 

5.250 E-01 

8.55 E-01 

Airframe 

-1.757 E+00 

4.174 B+00 

4.529 E+00 

3.880 E-01 

3.95 E-01 

Airframe 

-1.757 E+00 

-4.174 E+00 

4.529 E+00 

3.880 E-01 

3.95 E-01 

Engine 

-7.259 E+00 

0. 



9.55 E-02 

Inlet 

-2.725 E+01 

0. 



2.54 E-02 

Inlet 

-3.125 E+01 

0. 



2.22 E-02 

Inlet 

-3.139 E+01 

0. 



2.21 E-02 

Closed-loop eigenvalues for coupled system 

Airframe 

-5.329 E-02 

0. 



1.30 E+01 

Airframe 

-7.954 E-01 

7.665 E-01 

1.105 E+00 

7.201 E-01 

8.71 E-01 

Airframe 

-7.954 E-01 

-7.665 E-01 

1.105 E+00 

7.201 E-01 

8.71 E-01 

Engine 

-1.540 E+00 

0. 



4.50 E-01 

Engine 

-1.213 E+00 

1.135 E+00 

1.661 E+00 

7.302 E-01 

5.71 E-01 

Engine 

-1.213 E+00 

-1.135 E+00 

1.661 E+00 

7.302 E-01 

5.71 E-01 

Airframe 

-3.430 E+00 

4.169 E+00 

5.399 E+00 

6.353 E-01 

2.02 E-01 

Airframe 

-3.430 E+00. 

-4.169 E+00 

5.399 E+00 

6.353 E-01 

2.02 E-01 

Engine 

-7.212 E+00 

0. 



9.61 E-02 

Inlet 

-2.693 E+01 

0. 



2.57 E-02 

Inlet 

-3.125 E+01 

0. 



2.22 E-02 

Inlet 

-3.135 E+01 

0. 




2.21 E-02 




TABLE 9.- SIGNIFICANCE MATRICES FOR COUPLED, CLOSED-LOOP SYSTEM 


[Matrix order is similar to table 5] 


Eigenvalue 1 


1.615E-01 

1.687E-04 

7.695E-C5 

1.021E-02 

1.121E-01 

8.990E-06 

7.233E-01 

2.706E-04 

4.237E-C6 

9,720E-03 

1.195E-01 

5.302E-05 

9.168E-C2 

2.744E-05 

4.301E-C8 

l,02t;E-04 

5.858E-03 

2.092E-05 

3.053E-01 

3.408E-04 

1.685E-07 

6.972E-04 

3.493E-02 

3.032E-04 

2.800E-09 

1,049E-03 

9.418E-t4 

2.712E-02 

2.379E-02 

6.933E-03 

5,782E-05 

9.823E-03 

6.775E-C3 

1.711E-02 

1.496E-01 

9.478E-04 

1.6O1E-05 

5.826E-03 

6.198E-03 

l,183E-02 

1.294E-01 

5.924E-03 

5.5P5E-0<* 

9.320E-03 

1.759E-C2 

2.855E-02 

3.358E-01 

4.589E-02 

Eigenvalue 2 





1.991E-01 

6.471E-01 

1.659E-01 

1.C62E+00 

8.389E-01 

3.532E-02 

8.830E-03 

1.028E-02 

9,045E-C5 

l.OOlE-02 

6.853E-03 

2,062E-03 

4.256E-02 

3.964E-02 

3.492E-05 

3.995E-03 

1.651E-02 

3.095E-02 

6.470E-02 

2.248E-01 

6.243E-C5 

1.246E-02 

4.493E-02 

2,047E-01 

l.lOOE-01 

6.489E-04 

1.542E-D2 

6.418E-02 

5.196E-02 

2.284E-03 

2.249E-03 

6.015E-05 

1.098E-03 

4.009E-04 

3.236E-03 

3.091E-06 

2.«502E-02 

1.357E-03 

3.821E-02 

1.055E-02 

1,064E-01 

7,347E-04 

3.771E-01 

9*907E-04 

4.951E-C2 

1.162E-02 

1.261E-01 

2.598E-03 

Eigenvalue 3 





1.991E-01 

6.471E-01 

1.659E-C1 

1.062E400 

e.389E-01 

3.532E-02 

B,830E-03 

1.02PE-02 

9.045F-05 

l.OOlE-02 

8.853E-03 

2.062E-03 

4.256E-C2 

3.964E-02 

3.492E-C5 

3.995E-03 

1.651E-02 

3.095E-02 

6.470E-02 

2.24PE-01 

6.243E-05 

1.246E-02 

4.493E-02 

2.047E-01 

1,100E-01 

6.489E-04 

1.542E-02 

6.418E-02 

5.196E-02 

2,284E-03 

2.249E-03 

6.015E-05 

1.093E-t3 

4.009E-04 

3.236E-03 

3.091E-06 

2.502E-02 

1.357E-03 

3.P21E-02 

1.055E-02 

1.064E-01 

7.347E-04 

3.771E-01 

9.907E-C4 

4.951E-C2 

1.162E-02 

1.261E-01 

2.598E-03 

Eigenvalue 4 





1.500E-02 

8.242E-02 

3.774E-C3 

1 .733E-02 

f .497E-02 

4,612E-03 

1.956E-03 

3.850E-03 

6.052E-C6 

4.804E-04 

1.7C6E-03 

7,92tE-04 

3.951E-02 

6.221E-02 

9.788E-06 

8.034E-04 

1.333E-02 

4,980E-02 

7.146E-03 

4,19RE-C2 

2.083E-06 

2.983E-04 

4.317E-U3 

3.921E-02 

l,A36E-02 

4.500E-04 

2.455E-02 

2.401E-03 

9.670E-03 

2,222E-03 

8.637E-04 

1.227E-04 

f .142E-t 3 

4.410E-05 

1.771E-03 

8.846E-06 

4.026E-02 

1.159E-02 

7.495E-01 

4.860E-03 

2.439E-01 

8.e09E-C3 

7.222E-0? 

1.007E-03 

1.156E-C1 

6.370E-04 

3.441E-02 

3.707E-03 
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TABLE 9«- Continued 


Eigenvalue 5 


6.190E-0I? 

2.764E-01 

8,506E-t2 

3.619E-01 

1.594E-01 

1.529E-02 

1.161E-02 

1.857E-02 

1.961E-04 

l,443E-02 

7.114E-03 

3.776E-03 

2.636E-02 

3.373E-02 

3.566E-05 

2.712E-03 

6.249E-03 

2.669E-02 

4.593E-02 

2,193E-01 

7.310E-t5 

9.701E-03 

1.950E-02 

2.024E-01 

^•762E-02 

4.182E-04 

1.276F-02 

e.l07F-02 

3.535E-02 

2.430E-03 

4.U8E-03 

1.639E-04 

3.843E-t3 

2.141E-03 

9.310E-03 

1.391E-05 

2.158E-02 

1.742E-03 

6.298E-02 

2.654E-02 

1.442E-01 

1.557E-03 

3.72PE-01 

1.45BE-03 

9.355E-02 

3.351E-02 

1.959E-01 

6.313E-03 

Eigenvalue 6 





6.190E-02 

2.764E-01 

8.506E-C2 

3.619E-01 

1.594E-01 

1.529E-02 

1.161E-02 

1.857E-02 

1.961E-04 

1.443E-02 

7.114E-03 

3.776E-03 

2.636E-02 

3.373E-02 

3.566E-05 

2.712E-03 

6.249E-03 

2.669E-02 

4.593E-02 

2.193E-01 

7.310E-05 

9.701E-03 

l*950E-02 

2.024E-01 

A.762E-02 

4.102E-O4 

1.276E-P2 

P.107E-02 

3.535E-02 

2.430E-03 

A.lieE-03 

1.639E-04 

3.843F-C3 

2.141E-C3 

9.310E-03 

1.391E-05 

2.15BE-02 

1.742E-03 

6.298E-02 

2,654E-02 

1,4426-01 

1.557E-03 

3.729E-01 

1.458E-03 

9.355E-C2 

3.351E-02 

1.959E-01 

6.313E-03 

Eigenvalue 7 





4.740E-02 

5.418E-01 

6.878E-C1 

9.007E-01 

3.i99E-02 

3.202E-02 

1.029E-O3 

7,489E-03 

3.264E-04 

7.390E-03 

2,937E-04 

1.627E-03 

1,587E-C4 

5.197E-04 

2.267E-C6 

5.307E-05 

9.856E-06 

4.394E-04 

9.763E-04 

1.193E-02 

1.641E-05 

6.704E-04 

1.086E-04 

1.177E-02 

9.973E-02 

6.965E-04 

3.068E-C>3 

1.200E-02 

4.589E-03 

1.543E-03 

1.774E-03 

5#617E-05 

1.902E-04 

6.523E-C5 

2.487E-04 

1.817E-06 

3,552E-04 

2.280E-05 

1.190E-04 

3.088E-05 

1.471E-04 

7.771E-06 

2.168E-02 

6.741E-05 

I6.244E-04 

1.377E-04 

7.059E-C4 

1.113E-04 

Eigenvalue 8 





4.740E-02 

9.418E-01 

6.878E-D1 

9,007E-01 

3.199E-02 

3.202E-02 

1.829E-03 

7.489E-03 

3.264E-C4 

7.390E-03 

2.937E-Ci4 

1.627E-03 

1.587E-04 

5.197E-04 

2.267E-06 

5.307E-05 

9,856E-06 

A.394E-04 

9.763E-04 

1.193E-02 

1.641E-C5 

6.704E-C4 

1.086E-04 

1.177E-02 

9.973E-02 

6.965E-04 

3.068E-03 

l.POOE-02 

4.589E-03 

1.543E-03 

1.774E-03 

5.617E-05 

1.902E-C4 

6.523E-09 

2.487E-04 

1.817E-06 

3.552E-04 

2.280E-05 

1.190E-04 

3.088E-05 

1.471E-04 

7.771E-06 

2,16BE-02 

6.741E-05 

6.244E-04 

1.377E-04 

7.059E-04 

l,113E-04 
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TABLE 9.- Concluded 


Eigenvalue 9 


3.706E-04 
5.942E-03 
7.091 E-05 
3.945E-05 


9,3e6E-0A 

5.392E-03 

5.147E-05 

1.06RE-04 


?.0??E-04 

1.495E-03 

4.117E-09 

2.271E-04 


2.0O0E-05 

1.2e2E-05 

1.476E-09 

1.463E-03 


3.7P1E-13 

3.9CflE-15 

4.799S-15 

3.589E-13 


3,3O3F-03 
2.257E-05 
* 5 , 540E-06 
7.7A2E-04 


3.349E-05 

1.123E-03 

6.268E-05 

1.675E-05 


3.321E-04 

9.235E-04 

2.155E-03 

1.035E-C2 


2.609E-14 

1.222E-15 

&.043E-15 

1.102E-14 


2.001E-04 

6.034E-06 

3.003E-06 

2.033E-05 


1.541E-C3 
3.038E-04 
2 .903E-07 
1 .900E-t7 


6.733E-05 

1.734E-t3 

1.494E-04 

7.082E-C5 


1.410E-C4 

2.306E-06 

9.211E-t-8 

1.083E-C6 

1.611E-65 

3.443E-05 

1.239E-C4 

1.056E-03 


1.235E-14 

3.405E-18 

1.451E-19 

1.287E-18 

1.214E-15 

4.37JE-17 

1.678E-16 

1.079E-15 


4.391E-C4 

7.789E-C8 

6.635E-10 

l.lOOE-08 

1.1036-05 

I2.557E-07 

;i.963E-07 

I2.358E-06 


1.510E-03 

5.15CE-03 

5.0R8E-06 

5.810E-06 

7.237E-04 

1.635E-03 

1.C64E-04 

4.291E-05 


3.7C0E-05 

1.047E-C5 

4.322E-07 

8.867E-C6 

2.406E-05 

4.51CE-06 

1.227E-05 

8.889E-05 


2.795E-15 

1.332E-17 

5.868E-19 

9.082E-18 

1.103E-14 

3.483E-17 

l.OllE-16 

5.526E-16 


9.903E-05 

3.037E-07 

2.675E-09 

7.737E-08 

5.163E-05 

1.049E-07 

6.090E-06 

6.219E-C7 


2.364E-05 

9.020E-05 

4.164E-07 

4.148E-07 

4.598E-05 

1.035E-03 

e.425E-C5 

3.655E-05 


e.349E-oe 

2.643E-08 

5.100E-09 

9.126E-08 

3.6786-06 

6.873E-06 

2.338E-C5 

1.822E-04 


5.404E-17 

2.883E-19 

5.934E-20 

8.011E-19 

1.171E-14 

3.687E-16 

1.338E-15 

7.867E-15 


4.152E-07 

1.425E-09 

5.865E-11 

1.480E-09 

6.107E-05 

1.237E-06 

8.975E-07 

9.862E-06 


6.492E-05 

1.371E-03 

5.093E-05 

1.233E-C4 

2.363E-03 

1.157E-03 

6.806E-04 

8.809E-04 


6.712E-06 

1.176E-C5 

1.626E-05 

7.943E-04 

1.083E-05 

4.397E-07 

1.081E-05 

2.514E-04 


5.949E-14 

1.756E-15 

2.909E-15 

9.546E-14 

2.216E-15 

1.516E-18 

3.978E-17 

6.979E-16 


1.089E-03 

2.069E-05 

6.854E-06 

4.203E-04 

8.750E-06 

3.853E-09 

2.021E-08 

6.626E-07 


Eigenvalue 10 


1.467E-07 

1.952E-07 

9.735E-08 

9.729E-07 


1.743E-05 

8.305E-06 

3.314E-06 

1.236E-04 


Eigenvalue 11 

2.46PE-iy 7.312E-28 
5.531E-17 5.871E-30 
2.943E-17 2.499E-30 
2.219E-16 7.030E-29 


Eigenvalue 12 


1.996E-05 

2.R79E-07 

3.062E-08 

4.315E-07 


6.313E-05 

3.262E-07 

2.776E-C8 

1.459E-06 
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TABLE 10.- CLOSED-LOOP EIGENVALUES FOR DIFFERENT P 


Re(\) 

Im(X) 

Re( X) 

Im( X) 

Re(X) 

Im( X) 

(a) p 

= 0. 

(b) p 

= 0.01 

(c) p 

= 0.1 

-5.329 E-02 

0. 

-5.302 E-02 

0. 

-4.385 E-02 

0. 

-7.954 E-01 

7.665 E-01 

-7.987 E-01 

7.623 E-01 

-9.374 E-01 

6.938 E-01 

-7.954 E-01 

-7.665 E-01 

-7.987 E-01 

-7.623 E-01 

-9.374 E-01 

-6.938 E-01 

-1.540 E+00 

0. 

-1.532 E+00 

0. 

-1.507 E+00 

0. 

-1.213 E+00 

1.135 E+00 

-1.208 E+00 

1.139 E+00 

-1.114 E+00 

1.207 E+00 

-1.213 E+00 

-1.135 E+00 

-1.208 E+00 

-1.139- E+00 

-1.114 E+00 

-1.207 E+00 

-3.430 E+00 

4.169 E+00 

-3.431 E+00 

4.166 E+00 

-3.460 E+00 

3.993 E+00 

-3.430 E+00 

-4.169 E+00 

-3.431 E+00 

-4.166 E+00 

-3.460 E+00 

-3.993 E+00 

-7.212 E+00 

0. 

-7.197 E+00 

0. 

-7.230 E+00 

0. 

-2.693 E+01 

0. 

-2.692 E+01 

0. 

-2.726 E+01 

0. 

-3.125 E+01 

0. 

-3.125 E+01 

0. 

-3.125 E+01 

0. 

-3.135 E+01 

0. 

-3.134 E+01 

0. 

-3.131 E+01 

0. 

(d) p 

= 0.5 

(e) p 

= 1.0 

(f) p 

j = 00 

-5.123 E-02 

0. 

-2.191 E-03 

0. 

2.941 E-02 

4.471 E-02 

-6.418 E-01 

6.968 E-01 

-3.708 E-02 

0. 

2.941 E-02 

-4.471 E-02 

-6.418 E-01 

-6.968 E-01 

-4.465 E-01 

0. 

-9.664 E-02 

0. 

-1.394 E+00 

0. 

-9.264 E-01 

0. 

-4.443 E-01 

0. 

-8.750 E-01 

1.308 E+00 

-8.106 E-01 

1.314 E+00 

-8.106 E-01 

1.314 E+00 

-8.750 E-01 

-1.308 E+00 

-8.106 E-01 

-1.314 E+00 

-8.106 E-01 

-1.314 E+00 

-3.719 E+00 

4.599 E+00 

-7.259 E+00 

0. 

-1.757 E+00 

4.174 E+00 

-3.719 E+00 

-4.599 E+00 

-1.294 E+00 

7.374 E+00 

-1.757 E+00 

-4.174 E+00 

-7.212 E+00 

0. 

-1 .294 E+00 

-7.374 E+00 

-7.259 E+00 

0. 

-2.725 E+01 

0. 

-2.725 E+01 

0. 

-2.725 E+01 

0. 

-3.125 E+01 

0. 

-3.125 E+01 

0. 

-3.125 E+01 

0. 

-3.141 E+01 

0. 

-3.138 E+01 

0. 

-3.139 E+01 

0. 




TABLE 11.- CLOSED-LOOP EIGENVALUES WITH GAIN ELEMENT (1,3) REMOVED 


Subsystem 

Re(\) 

lm( X) 

03 

c 

^1/2 ^2 

Airframe 

-5.328 E-02 

0. 



1.301 E+01 

Airframe 

-9.024 E-01 

6.250 E-01 

1.098 E+00 

8.221 E-01 

7.681 E-01 

Airframe 

-9.024 E-01 

-6.250 E-01 

1.098 E+00 

8.221 E-01 

7.681 E-01 

Engine 

-1.535 E+00 

0. 



4.516 E-01 

Engine 

-1.113 E+00 

1.149 E+00 

1.600 E+00 

6.958 E-01 

6.228 E-01 

Engine 

-1.113 E+00 

-1.149 E+00 

1.600 E+00 

6.958 E-01 

6.228 E-01 

Airframe 

-7.404 E-01 

5.607 E+00 

5.656 E+00 

1.309 E-01 

9.362 E-01 

Airframe 

-7.404 E-01 

-5.607 E+00 

5.656 E+00 

1.309 E-01 

9.362 E-01 

Engine 

-7.207 E+00 

0. 



9.618 E-02 

Inlet 

-2.694 E+01 

0. 



2.573 E-02 

Inlet 

-3.125 E+01 

0. 



2.218 E-02 

Inlet 

-3.134 E+01 

0. 



2.212 E-02 
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